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Using a common technique for approximating distributions [generalized functions], we are able to use 
standard Monte Carlo methods to compute QFT quantities in Minkowski spacetime, under phase transitions, or 
when dealing with coalescing stationary points. 
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1. Theory 

We use the mollification {approximate identity) 
technique to rewrite the partition function in such 
a way that not only its computation, but the com- 
putation of its derivatives (Green's functions), be- 
come well defined using Monte Carlo methods, 
even for integrands that normally do not allow 
this approach. 

1.1. Mollification 

A mollifier, rj, is a positive, C°°[/] function, 
(/ G M), with supp(77) = B{OJ), where I = 
length(/) and B{0, 1) is the ball centered at 
with radius / and Jjil = 1 • The sequence of func- 
tions, rye = e~"ri{x/e) is an approximate iden- 
tity. A mollification of a locally integrable func- 
tion / : — > M, ([/ e M", open), is given by. 



Ue = {xeU\ dist(a;, dU) > e} 
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B(0,e) 



■q^iy) f{x -y)dy . 



1.2. Properties 

The properties of mohifiers, (/e), are: 

1. /e eC°°(C/,); 

2. fe ^ f, almost everywhere, as e 0; 



3. If / £ C{U), then fe f uniformly on 
compact subsets o{U\ & 

4. If 1 < p < cx) and / e L\^^{U), then /, -> / 

It is not difficult to see that this approximation, 
(/e, the mollified version of /), can smooth oscil- 
lating functions, /c is differentiable and, not only 
does it converge to / when e ^ but, depending 
on where this convergence takes place, it can be 
either uniform or in L^^^- 

1.3. Implementation of the Method 

The object we calculate will not be the usual 
partition function, but its mollified version, Ze, 



where. 
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(/), s /, ^ 77, * / . 

From equation (j|) it follows that, upon choosing 
U appropriately, the boundary condition^ for the 
theory (defined here by S) are set|^. 
The method is apphed in two steps: 
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^Note that, a choice of boundary conditions [Schwinger- 
Dyson equation] is equivalent to choosing the measure in 
the path-integral. 
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1. Saddle-point expansion : For simplicity, the 
mollifier is chosen to be a Gaussian^, 

^-(^^ ^ V(2^)" dct(e^) ' 

where is the covariance matrix that de- 
fines the Gaussian distribution. Thus, a 
typical path-integral will look like, 

3e = j {f{v>)e'-'^'^^-''-'^)^ m 

oxpj -| -{l + e^ ■ S" ■ 6)-i ■ b| 

X f -r ^= M + 

-f corrections , 

where B = e • VS{ip) and = 
d'^S{ip)/dipidipj (Hessian matrix for S). 

2. Importance sampling : The Monte Carlo 
simulation is handled as follows: Choose an 
Importance Function, {W{ip)), and change 
the measure of integration accordingly, 



J 



i S{(p) — i J-ip 



J 



[dW{ip)] 



(5) 



W,(¥>) 



A reasonable choice for W is given by | 
W{ip) = Weiif) = Ke'^^"^^),!- Upon 
saddle-point expansion, we find, 

[iSiifio) - I B"r ■ (1 -h e"r ■ S" ■ e)-^ ■ b} 



cxp< 



ydct(l + eT ■ S" ■ e) 



(6) 



2. Applications in Physics 



Two basic examples will be shown here. Both 
of them can be regarded as 0-dimensional [space- 
time] quantum field theories. 

^Note that, bold symbols and bold letters, denote vec- 
tor/matrix quantities. 



2.1. Airy Action 

The action is given by: S{x) — ix^/3 + itx. 
Note that, exp{S'(a;)} dx = Ai(t), thus, the 
partition function and its mollification are. 
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3 ^ Ai(0) 
dx 
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+ itx>) dx 



(8) 



I — > } dx 



This is a particular useful example for two rea- 
sons: 

1. We can easily compare the approximated 
results with the exact ones, 



2. Boundary Conditions: Assuming integra- 
tion on the real axis, the Euclidian version 
of this theory does not exist, therefore work- 
ing in Minkowski space is necessary. 

It is easy to see that this is a highly oscillatory 
integrand (what makes the use of Monte Carlo 
techniques impractible). However, once its mol- 
lification is calculated, the task becomes simpler 
(as shown below). 




! nihil 



The importance function, (W(x)), is also shown 
for several values of the parameter e: 
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-8-6-4 -2 2468 
Legend 



W(x,-1 6,0.30) 
W x,-16,0.08 
Wx,-1 6,0.38 
W x,-16,0.68 



2.2. cjy- Theory 

The action is, S{4>) = ficj)'^ /2 + g (jj^ /A. Thus, 
the molUfied generating functional is given by: 




As usual, its parameters, {rn,g), will determine 
the phase structure of the theory. The follow- 
ing are the plots of the 2-point Green's function 
in the symmetric and broken-symmetric phases, 
respectively. 




The plots below are the importance function and 
the mollified integrand in the symmetric and bro- 
ken symmetric phase, respectively. 




It is interesting to reahze that, for the broken 
symmetric phase, the importance function peaks 
exactly at the point where the mollified integrand 
ceases to have two concavities and begins to have 
only one. 
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